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Abstract 

A general formalism for the study of excitations above equilibrium in Bethe ansatz 
solvable models is presented. Nonzero temperature expressions for dressed energy, mo- 
mentum, spin and charge are obtained, and it is found that the dressed spin and charge 
are in general momentum dependent. For an electronic model one may only have spin- 
charge separation at zero temperature where the ground state is half-filled and has zero 
magnetisation. Finally, the excitations of the Hubbard- Shastry models are examined in 
detail. 
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1 Introduction 

The quest to understand strongly correlated electrons in low-dimensional systems repre- 
sents an important frontier in the field of condensed matter physics [HE]. Such materials 
exhibit exotic behaviour that cannot be understood from a non-interacting picture. Ex- 
actly solvable models are an invaluable tool for gaining access to the non-perturbative 
physics at play. 

Two recently introduced integrable models, the Hubbard-Shastry A- and B-models 
[3], hold great promise in this direction. They describe electrons interacting on a lattice, 
similar to the well known Hubbard model [U E] , but with extra interactions such as spin 
exchange, pair hopping and nearest neighbour Coulomb interaction. Their equilibrium 
properties were examined in [3] by means of their exact Bethe ansatz solution. It was 
found that the A-model exhibits itinerant ferromagnetism, while the B-model is a Mott 
insulator of paired electrons and so provides a promising approach to the study of un- 
conventional superconductivity given that such physics is expected to be captured by 
an effective single-band model [HI E] • 

The aim of this paper is to examine the excitations above equilibrium of these two 
models in one- dimension. To achieve this goal a general picture of excitations in Bethe 
ansatz solvable models is pursued. In particular our consideration leads to natural 
expressions for dressed energy, momentum, charge and spin, and from these some general 
conclusions can be made. Most noteworthy is that for electronic models, including the 
Hubbard model, the spin and charge of excitations above equilibrium get dressed and in 
general depend on the momentum of the excited quasi-particles. The spin and charge 
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scale with the energy of the excited quasi-particle, see for example figures [9 10 This is 
possible as excitations also give rise to an induced spin and charge of the system that 
does not propagate. The spin and charge only become constant at high temperature 
where they take their undressed values, and in certain regimes at zero temperature that 
correspond to half-filling and zero magnetisation. Away from half-filling, spin and charge 
are in general not separated and so one may expect a breakdown of the Luttinger liquid 
picture, see e.g. [8]. 

These features are a consequence of the strongly interacting nature of such models. 
We would like to suggest that non-trivial dressing of charge and spin may account for 
some of the anomalous behaviour observed in strongly correlated materials. It would be 
interesting to consider issues such as the unusual temperature dependence of resistivity 
and Hall coefficient in this context, see e.g. [9]. 

The paper is naturally split into two parts. The first provides a general formalism 
for the study of excitations above equilibrium in Bethe ansatz solvable models. This 
generalises and extends the methods used in [TUJ-JTS] . In particular the dressing of 
charges corresponding to conserved quantities is considered. In the second part of the 
paper the formalism is used to examine zero temperature excitations of the Hubbard- 
Shastry A- and B-models. 

Bethe ansatz solvable models represent an important class of one-dimensional models. 
In them the complexity of calculating the spectrum is reduced to the solution of the 
Bethe equations, a set of polynomials whose degree scales linearly with the length of the 
system. In the thermodynamic limit the roots of these equations, in general complex, 
align into strings on the complex plane known as Bethe strings. It is often possible to 
make a string hypothesis, which identifies all possible Bethe strings, and in this paper 
we restrict our attention to models for which this is so. Corrections to the Bethe strings 
are exponentially suppressed in the length of the system. 

There are advantages to working directly with the Bethe strings in the thermody- 
namic limit. The Bethe strings are parametrised by a rapidity variable that is real 
whereas, as we have mentioned, their constituents (the roots of the Bethe equations) 
generally take complex values. The original Bethe equations can be fused together to 
give Bethe equations directly for the strings. Moreover, the Bethe strings give a physical 
picture: the different types of strings can be understood as the particle content of the 
model. The free energy can thus be written and its minimisation gives access to the 
equilibrium state. In this way one goes from the full information about the spectrum 
that the Bethe equations provide to the physically interesting region. 

To examine excitations above equilibrium it is also natural to work directly with the 
Bethe strings and this is the approach taken in the formalism developed in section [2] 
Excitations of a finite number of Bethe strings are considered. One finds that when 
excited above the equilibrium state the energy, charge and spin of the Bethe strings get 
dressed. It is noteworthy that the dressed charge and spin are, in general, rapidity de- 
pendent. Moreover, the scattering of excited roots is non-trivial and includes scattering 



with the equilibrium state. A summary of these results is provided in section 2.5 



In the zero temperature limit the equilibrium state becomes the ground state. This 
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limit is in general smooth, with the special exception of when some strings are at half- 
filling. By a half-filled string we mean one for which there are no holes in the ground 
state. In such cases there may exist restrictions on allowed excitations over the ground 
state. This can be seen as a manifestation of a zero temperature phase transition. 

Now we turn to our study of the excitations of the Hubbard-Shastry models, which 
comprises the second part of this paper. Attention is focused on the zero-temperature 
regime as the essential features can be observed there. General features of finite tem- 
perature behaviour are discussed throughout the presentation of section [2} A detailed 
analysis of the finite temperature excitations requires the solution of various infinite sets 
of coupled non-linear integral equations and is not pursued here. 

The A-model has a ground state which is spin polarised. At half-filling there are 
magnons which disperse quadratically and so one does not expect spin-charge separation. 
There is an electron-like excitation which is gapless away from half-filling and becomes 
gapped at half-filling. The excitations that behave as magnons at half-filling retain their 
quadratic dispersion away from half-filling but here get dressed with some charge in 
addition to their spin. The dressed magnons are gapped in the presence of a magnetic 
field. 

The B-model has a richer ground state phase diagram. At half-filling the model is in 
an insulating phase. If there is no magnetic field then the magnetisation is zero and the 
model exhibits spin-charge separation. All excitations are scattering states of gapless 
spinons and gapped holons and the picture is similar to that of the Hubbard model [3]. 
The difference is that here the spinons are dressed electrons and the holons are dressed 
spin singlets (paired electrons) whereas in the Hubbard model the spinons are dressed 
magnons and it is the holons that are dressed electrons. If one goes away from half- 
filling or introduces a magnetic field then the spin and charge are no longer separated 
and become rapidity dependent. In particular away from half-filling in the absence of 
a magnetic field the dressed spin singlets remain pure charge and become gapless while 
the electrons are also gapless but get dressed with a charge in addition to their spin. 
Thus one may expect behaviour beyond that of a Luttinger liquid. The dressed charge 
of an excited electron does go to zero however as its dressed energy goes to zero. Let us 
conclude by remarking that the dressed electrons display an "hourglass" dispersion and 
that away from half-filling and in the absence of a magnetic field it can be clearly seen 
that the spin singlets are held together by spin-spin interactions. 



2 Excited states 

In this section we develop a general formalism to study excited states of any Bethe ansatz 
solvable model whose thermodynamics is described by various configurations of Bethe 
strings [H]. The dressed energy, momentum and charge, and also the phase shifts of the 
excitations, are expressed through the densities of the equilibrium state. The limit of 



zero temperature is examined in section 2.4 and the formalism is summarised in section 



2.5 For simplicity and clarity of presentation we restrict ourselves to parity invariant 
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rational models. 



2.1 Bethe equations for strings 

Let us assume that in the thermodynamic limit every solution of the Bethe equations 
for the model under consideration corresponds to a particular configuration of Bethe 
strings. Then for large spin chain length L the Bethe equations can be written for string 
configurations 

= e iL Pa (v a , k ) "Q jj S a0 (v a)k ,v p>n ) . (2.1) 

P n=l 

We use indices a, (3, 7 and so on to distinguish between different strings. Here ip a is a 
constant, unimportant for our consideration, which appears in particular as we have not 
excluded self-scattering from the product on the right hand side. For periodic boundary 
conditions and in the absence of a twist, (p a = if S aa (0) = 1 and <p a — 1 if S aa (0) = — 1. 
The strings are parametrised by a rapidity variable for which the scattering matrices are 
of a difference form: S a p(v at k,Vp )n ) = S a p{y a) k — vp,n)- For strings which do not carry 
momentum there is no term e lLpa in the Bethe equations, and by convention we set 
p a = for such strings. 

Let us briefly review the analysis that yields the equilibrium properties of the model. 
This allows us to introduce many useful formulae, along with our conventions. The 



counting function is constructed by taking the logarithm of the Bethe equations (2.1) 
La a z a (v) = irip a + Lp a (v) + ^2^2 </WK u /3,n) , (2.2) 

P n=l 

where ^ 

<t>ap = T log Sap (2.3) 

is the phase of the S-matrix. These functions allow one to enumerate the Bethe strings 
as h^- evaluates to an integer on an a-string Bethe root, which we refer to as the 
mode number of the particle. For a given root v a ^ we denote the corresponding mode 
number as I a ,k = j^z a {v a ^)- Moreover the counting function Lz a /2n may evaluate to 
an integer for a momentum which is not that of a particle of type a, and such mode 
numbers correspond to holes. Note the appearance of a a in the definition of the counting 
function. For momentum carrying strings a a = sign(^) and this guarantees that the 
counting function is an increasing function of v. For an a-string which does not carry 
momentum a a is determined by requiring the counting function to be increasing. Let 



us remark that to fully define the counting function (2.2) it is necessary to specify that 
branch of the logarithm. This does not affect the study of equilibrium properties, and of 
course should not affect the physics, but it will be important for our study of excitations 



and we will return to this issue when we begin to discuss them in section |2.3 
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Taking the thermodynamic limit L — > oo with N a /L fixed, one gets equations for 
the densities of particles and holes 



Pa + Pc 



1 

2^ 



dp a 



dv 



+ K a f3 -k pfs , 



(2.4) 



where Lp a dv (respectively Lp a dv) is the number of integers corresponding to particles 
(respectively holes) that the counting function Lz a /2ir evaluates to over a range dv. 



Repeated indices are summed over and * denotes convolution (see appendix A.l for the 



precise definition) over the domain of the rapidity of the appropriate string, which we 
denote by I a . The kernels K Q/3 are defined by 



K, 



a/3 



a/3 



K a p{v) 



1 d 

2i\i dv 



log S a p(v) 



(2.5) 



Since the counting functions are all defined to be increasing functions it follows that the 



densities of all particles and holes in the equations (2.4) are positive 



The equations for the densities can be used to determine the total number of particles 
and holes for each type of string, which we call the range of mode numbers. Indeed, 



integrating equations (2.4) and multiplying by L one gets the range 



N a + N a 



L — h k a pNp 

2n 



(2.6) 



where N a is the number of holes of type a in the state, and Ap a = 1*| -J^ | for momentum 
carrying strings and Ap a = for strings that do not carry momentum. The k a @, defined 
as k a p = 1 * K a p, are constants. Note that all roots of the Bethe equations are counted 
in the range, including those which correspond to v = oo. Furthermore since N a > 



eqs. (2.6) imply the following selection rules 



N a < L ^ + KpNp , 
which restrict the allowed N a appearing in the Bethe equations. 



(2.7) 



2.2 Equilibrium state 

The equilibrium state follows from the minimisation of the free energy 

/ = e a * p a - Ts . 



(2.1 



Here e a denotes the energy carried by an a-string. Note that in addition to describing 
possible dispersion this will also depend on the chemical potentials that appear in the 
model. The entropy per site is denoted by s and is given by 



s = Y] s(p a , p a ) , s(p, p) = p log (l + -)+p log (l + -) • 
' P P 



(2.9) 
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Minimising the free energy subject to the equations for densities (2.4) yields 



/ 



T 
2k 



dv 



*log (1 + — ) 

J- rv 



(2.10) 



where the functions Y a are determined by the set of thermodynamic Bethe ansatz (TBA) 
equations 

log Y a = % - log (1 + ^) * Kp a . (2.11) 

Y I3 



These should be regarded as the conditions for equilibrium. The Y a are related to the 

Pa 



densities of particles and holes as Y a 



Pa 



and allow one to solve (2.4) for the densities 



at equilibrium. It is useful to also introduce the pseudo-energy 

e a = Tlog Y a . 



(2.12) 
Both 



These are well defined in the zero temperature limit whereas the Y a are not. 
quantities Y a and e a will be useful and we will use them interchangeably. 

At finite temperature the equilibrium state is a mixture of infinitely many states 
contributing equally to the partition function. We assume for definiteness that the spin 
chain length L is very large, and choose any of the states which in the thermodynamic 
limit contribute to the equilibrium state as a reference state. In the zero temperature 
limit it becomes the ground state of the model. Let us say then that reference state 
consists of N a a-strings with rapidities v a k, k — 1, . . . , N a . 



2.3 Excited states 

Now we wish to study excitations about the equilibrium state outlined above. We restrict 
our attention to excitations where the numbers of excited roots are much smaller than 
the numbers of particles in the equilibrium state. 

Let us first return to a point that we skipped above, the choice of the branch for the 



counting function (2.2). The choice of branch affects the formalism one obtains for the 
excitations and we present here what we found to be an optimal choice. In particular 
we find that in order to obtain a reasonable expression for the dressed momentum of 



a string it is necessary to keep track of the branch of each log S a p term in eq.(2.2). 
Moreover, our choice of branch is guided by the behaviour of the pseudo-energies for the 
equilibrium state. We choose the branch of log S a p and the range of momentum p a so 
that the counting function is continuous about the minimum of e a . 

Let us be specific about the strings we will consider. For strings with rapidity variable 
defined on K we assume that the pseudo-energy is even, and monotonic on the interval 
(0, oo). Then there are two cases, and we give explicit expressions for the phase in each 

• Type 1: e a (v) has a minimum at v = 0, and is increasing on the interval (0, oo), 

T log S a/3 (V, t) = 2ltb a p + 7TC a/3 + <d aP (v - t) = <j>ap(v, t) , (2.13) 

and the range of p a (v) is chosen so that it is continuous for v along (—00, 00). 
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Here 



Type 2: e a (v) has a minimum at v — ±00, and is decreasing on the interval (0, 00), 
-logS a/3 (v,t) = 2irb a/3 - 7rc Q/3 sign(f) + 9 a p{v-t) = </>ap(v,t) , (2.14) 

% 

and the range of p a (v) is chosen so that p Q (— 00) = p a (+oo), and it is discontinuous 
only at v — 0. 

f v 1 

Q a p{v) = 2tt / dtK, a p(t) , = -6^(00) = l*K a p , (2.15) 

is an integer capturing the freedom in the choice of branch, and the c a/ 3 are defined 
such that c a/ 3 — is as close to zero as possible subject to the constraint that c a/ 3 is 
an even integer if S a p(0) = 1 and c a p is an odd integer if S a p(0) = —1. Note that if 
S a p(±oo) = 1 (which is the case for rational S-matrices) then c a/ 3 = K a p. For type 1 
strings </> a ^(±oo,t) = 7r(c a/ 3 ± K a p) and the range of <fr a p is Qj g(+oo,i) — a( g(— 00, t) = 
27TK a /3. For type 2 strings it is worth stressing that the scattering phase Q/ g is no longer 
of a difference form with this choice of the branch of logS a p. It is mildly broken so 
that, for models with c a p = K a p, the counting function and scattering phases would be 
continuous everywhere but at v = 0, the maximum of e a . The jump discontinuity of the 
scattering phase at v = is equal to (j) a p(—0,t) — (j) a p(+0,t) = 27rc a/ 3, which is equal to 

27TK a/ g if S a p(±Oo) = 1. 

These two cases capture the behaviour of most strings of Bethe ansatz solvable 
models. For strings that are not captured modifications of the subsequent analysis will 
have to be made. In particular, the y-particles of the Hubbard-Shastry models must be 
considered separately because their rapidity variable is not defined on R. We deal with 



this important case in appendix |A.3 



Now consider a general excited state. Let us say that it consists of N a a-strings with 



rapidities v a k, k — 1 , . . . , N a . The rapidities satisfy the same Bethe equations (2.1 ) 



(-l)v- = e u*«(««, h ) Y[ J] S a p(v a>k - v p , n ) . (2.16) 

n=l 

The rapidities v a ^ of the excited state can be divided into two groups. The first group 
consists of rapidities with mode numbers which coincide with some of the mode 
numbers of the particles of the reference state. They are close to the corresponding 
rapidities of the reference state, that is, the difference between the rapidities with the 
same mode number is of order 1/L. We denote these rapidities as v' ak , k = 1 , . . . , N' a . 
The second group consists of the remaining rapidities, those which have mode numbers 
not coinciding with any mode number of the particles of the reference state. There are 
= N a — N' a of strings of type a in this group. These rapidities will be denoted as 
v a .. The reference state also contains strings with mode numbers different from any 
I a> k, the mode numbers of the particles of the excited state. These correspond to holes 
of the excited state and their rapidities will be denoted as v x .. For a string of type 
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a there are iV* = N a — N' a of them. One can think about the excited state as being 
obtained by adding N* strings to, and removing strings from, the reference state. 
Let iV a = J2 a be the total number of strings added to the reference state, and let 
N r = Y2 a be the total number of strings removed from the reference state. Thus the 
equations (2.16) can be rewritten in the form 

= e — k,j njj^y n n Sapi<k _ vkn) > (217) 

LLj ^atj[V aik V t .) p n=1 

where the product n^IIn^i Sap('> v 'p,n) includes the product Ylf £>atj(',v tj ), that is it 
is equal to fL Iln=i ^(-^kj IlT S av (-,v t ) and reduces to fLECi S a p(-,v 0>n 



in 



the thermodynamic limit. 

Here we have made an implicit assumption that the mode numbers of every string 
in the reference state are also mode numbers of the excited state. This is justified at 
non-zero temperature as one can always choose the reference state so that this is the 
case^ At zero temperature however, specifically at half-filling, it may happen that the 
mode number of a string in the reference state is not a mode number of the excited 
state, either due to a change in the range of mode numbers or an overall shift of the 
range of the counting function. For example if the range is decreased for an excitation 
and all mode numbers correspond to particles then some strings are necessarily removed, 
and moreover there is no excited state rapidity one can assign to them. Such situations 



require special care and are discussed in section 2.4 



We relate the rapidities of the reference state Bethe roots to those with corresponding 
mode number in the excited state through 

<3a Ca(Pot,k) / 1r .N 

V a ,k ~ V a)k = , (2.19) 

Z7T L 

where we have introduced ( a which are of order 1. It is possible to obtain a closed 
equation for ( a by subtracting the logarithm of the Bethe equations of the ground state 
(2.17) from those of the excited state (2.1) for Bethe roots with the same mode number. 
Expanding (2.17) and taking the thermodynamic limit one obtains 

(a{Pa + Pa) = CpPfi * Kp a ~ <f) aSu + <f) ar , (2.20) 



with the help of the equation for densities (2.4). Here we have taken into account that 

Sap{v)Sp a {-v) = 1 => ^^-]og3 a p{v-t) = -K Pa {t-v), (2.21) 

Lixi at 



1 Indeed the densities p(v), p(v) are related to the numbers of particles and holes with rapidity in an 
interval dv about v. At non-zero temperature both the densities of particles and holes are non-trivial 
for all v and the particles and holes within each interval dv can be rearranged. Thus one can always 
choose a reference state at non-zero temperature such that a finite number of particles and holes have 
required mode numbers. 
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and introduced the notation 



JV a N' 
3=1 3=1 



(2.22) 



for any quantity X a . The function F a = —( a (pa + Pa), which we will refer to as the shift 
function [HI [12], is an important object and it satisfies the following closed equation 



1 P 



*K dQ + 



0a T (Pan — H>ax 



u 1+Yp 

Another form of this equation which will prove useful is 



(2.23) 



K 



Energy 

The change in the energy of the excited state from the equilibrium state is given by 



■'aa. far 



(2.24) 



(2.25) 

/ Z7T 

a k=l k=l 

where e' a = ^-e a (v), and summation over a is assumed, and we have used the notation 
(2.22). Differentiating the TBA equations (2.11) and integrating by parts one gets 

6 f = (log YJ -IU (2.26) 



and substituting into (2.25) gives 



AE = e a - e r + T°^( aPa * ((log Y a ) f - * (log^)') 

Z7T 1 + Yr 



(2.27) 



Dependence of AE on £ Q can be eliminated through eq. (2.24) yielding 

T 1 
AE = e a - e r + — <r tt (log(l + — )) * (0 aa - ar ) . 

Then integrating by parts gives 

A£ = e a -e r -T log(l + — ) * (A aa - A Qr ) + Tlog(l + ^^){K, - k aI ) , 



(2.28) 



(2.29) 



2 Let us remark that the jump discontinuity of 4> a p for type 2 strings is 2irc a p, which is equal to 
2-Kkap only if 5 Q| g(±oo) = 1. 
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where Y^ nax = Y a (v max ) is equal to Y a evaluated at the value of v corresponding to the 
maximum of the pseudo-energy, and recall k a p = l-kK a p = (T a K a p- Finally we use again 



the TBA equations (2.11) to obtain 



7V a JV r 



3=1 $ ; i 

Hence the dressed energy of an a-string is 



(2.30) 



E a {v) = e a (v) + Tlog(l + Y^)k Pa 



(2.31) 



Note that the last term in the formula is rapidity independent and does not contribute 
to the total energy in a particle-hole excitation. Let us mention that in limit of infinite 
temperature the dressed energies take their bare values E a = e a as one should expect. 
This follows as the functions Y a become constant in the limit T — > oo because the driving 



terms drop out of the TBA equations (2.11) 



Momentum 

In a similar way the change in the momentum of the excited state from the equilibrium 
state is given by 



a k=l k=l 



dp 



dv 



(2.32) 



Let us first remark that the momentum of a state is defined modulo 2n and thus, as 
the momentum of the reference state is fixed, the change in momentum is also defined 



modulo 2tt. To simplify expression (|2.32|) we substitute (]2.4|) into (|2.20|) and get 



2tt 



(a 



dv 



aa i far 



(2.33) 



Multiplying by p a , integrating, and taking the sum over a we find 



2tt 



Cot Pa ~k 



dp c 



dv 



-Pa * 



'oa far J j 



and thus 



Hence we identify the dressed momentum of an added a-string as 



)• 



(2.34) 



(2.35) 



Pa = Pa ~ P/3 * 



(2.36) 
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and a removed one with the opposite sign. The dressed momentum can be used to 



parametrise the strings and this is discussed in appendix A. 2 



To examine the range of dressed momentum it is useful to note that 



dP a 



dv 



2n(p a + p c 



(2.37) 



which is seen using eqs. ( |2.4 ) and (2.21). For a-strings of type 1 the range of dressed 
momentum is over (P Q (0) — n(n a + n a ),P a (0) + Tc(n a + n a )) where P a (0) = p Q (0) + 
27inpbp a + irnpcp a , n a = 1 * p a and n a = 1 * p a . On the other hand for strings of 
type 2 the range is split into two parts. Recall that the bare momentum of a type 2 
a-string has a jump at v = and that it increases from p Q (cr Q 0) to p a (— & a 0). The range 
of dressed momentum in this case is thus over (P a (cr Q ,0), P a (a a 0) + n{n a + n a )) and 
(P Q (— (JqO) — vr(n a + n a ), P a (—a a 0)). In general this may result in a gap in the dressed 
momentum. 

Let us remark that the dressed momenta depend on the choice of branch of log S Q p, 
i.e. the b a p in eqs. (2.13), (2.14). Considering only strings of type 1 and 2 the dressed 



momenta (2.36) can be written as 



Po 



P/3 * - 2,Tm b 0o 



^2 7rn /3 C /3" > 
/3 of type 1 



(2.38) 



where one is free to choose the integers bp a . Terms for cp a with (3 of type 2 do not 
contribute as we restrict ourselves to parity invariant models and so the densities are 
even. It is of course possible to describe any excitation with a definite choice of bp Q , e.g. 
one can set all bp a = 0. However in this case some excitations would have unnatural 
description which would require considering particle-hole excitations with zero energy 
contributing only to the total momentum. For np irrational one can achieve any value 
of dressed momentum by choosing bp a appropriately. Let us remark that one is free to 
choose b/3 a independently for each added and removed a-string and one may refer to the 
set b/3 ad , bp Xj , where /3 runs over all strings that interact with the excited string, as the 
branch of the excitation. 



Phase shift 

Now we turn our attention to the scattering phase shift. Consider first the counting 
function for an a-string of the excited state 

L (J a Z a (v) = 7lip a + Lp a (v) + (fia^v) - (p aT (v) + ^ ^Z 3 ^' V P,n) ' ( 2 " 39 ) 

P n=l 

Expanding v'g n in the final term about its equilibrium value, replacing the sums by 
integrals, and noting equations (2.20, 2.36), one gets 



L cr a z a = TC(p a + LP a + F a . 



(2.40) 
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Recall that F a here is the shift function which is determined through the closed set 



of equations (2.23). Exponentiating equation (2.40) and evaluating it at a rapidity v 



corresponding to a mode number of the excited state it takes the form 



(2.41) 



An added a-string with rapidity v has dressed momentum P a (v) and so its scattering 
phase shift is 5 a = F a (v)+7T(p a . Similarly a removed a-string with rapidity v has dressed 
momentum — P a (v) and so its scattering phase shift is 5 a = —F a {y) — ir<p a - Clearly the 
phase shifts 5 are defined modulo 2n. 



As the equations (2.23) are linear it is natural to introduce the set of functions 
&ap(v,t) satisfying the following system of equations 



± r 



+ 1 + K 7Q ' 



(2.42) 



where it is understood that 



\ l+Y 1 - fv 7« 



[v,t) 



Then F n 



$ Qa — <3> OT and hence we refer to the $ a « as dressed scattering phases. In terms of these 



functions equation (2.41) takes the following physically intuitive form 



7V a N 1 
e iP a (v)L TT g»a j (v ,v aj ) JT e - 



(2.43) 



It is worth mentioning that the dressed scattering phases are in general not of a difference 
form, and in particular $ ai s(i>, v) ^ 0. This is a reflection of the fact that an excitation 
has nontrivial scattering with the equilibrium state. 



Charge 

Each chemical potential appearing in the model is related to a conserved quantity and a 
corresponding charge. For a given chemical potential fi let us denote the corresponding 
bare charge 3 carried by an a-string as w Q = —^g- Furthermore let us introduce an 
object uj a = —^7, which we call the pseudo-charge of an a-string. It satisfies the 
following set of integral equations 

w« = w a + * Kj3 a . (2.44) 

J- + ip 

The change in the total charge of the excited state from the equilibrium state is 

N a N a 

AW = ^ f y] ™ a (v a)k ) - ^ W «(^a,fc)) ^ W a - W r + ^(aPa * , (2.45) 
a k=l k=l 



3 For models where a magnetic field B enters as a chemical potential the bare spin of an a-string is 
generally given by 
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where here we are being formal as w Q has no rapidity dependence. Indeed the final term 
is zero, but let us further analyse it nevertheless. Recalling that ( a p a = —jtt^, we have 

J- i X a 

= -TTV a * UJ " + F «*(iTvJ ~^^1TVJ + ^^YTVJ 

1 / Ul y max 

=F a u a -k ( ) + 2ixa a a -k (K a3i - K ar ) - 2ncr a a max (fc aa - k aT ) . 

1 ~r I a J- ~r I a 1 ~r '„ 

Here cu™ ax = u a (v ma,x ) is defined similarly to Y™ a *, both functions being evaluated at 
the value of v corresponding to the maximum of the pseudo-energy. This allows one to 
write the total change in charge as 

J " . .max Jv . .max i 

AW = V (co a .(v a .) ^ k aa .)-Y (wt-fe.) — JW)+— W)'- 

/ j V a J V n J ' J_ _|_ ymax Qa J / / / V r J V r J / ^ _|_ ymax ax J I ^ a a \ ^ _|_ Y ' 

3=1 a 3=1 a ~ a 

We would like to present an interpretation of this change as 

iV a N T 

aw = j2 w aj ) - Yl w * ) + AWind • ( 2 - 46 ) 

i=i i=i 
Here the excited strings are assigned a dressed charge 

W a (v) = u a (v) - Y^^kpa = , (2.47) 



that they carry, while the final term 



In 1 + 



kh-y (2 - 48) 



is understood as an induced charge of the system that cannot be transferred. An in- 
teresting feature here is that the dressed charge carried by an excited string depends 
in general on the string's rapidity. In the limit of infinite temperature the functions Y a 
become constant and the dressed charges take their bare values while the induced charge 
goes to zero. The zero temperature limit will be discussed in the next section where it 
is seen that the induced charge resides at the edge of the Fermi sea. 

The above is the interpretation we shall adopt in this paper but let us mention that 
the final term in eq. (2.46) can be redistributed among the added and removed roots 
using F a = $ aa — $Q, r . In particular, integrating by parts this final term one obtains 
back 

AW = J2w aj -J2w tj , (2.49) 

3=1 3=1 
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via the curious identity 



1 + Y B 



*K, 



2tt1 + Y 6 



-k (Up$p a )' 



(2.50) 



One may wonder why we insist on the interpretation of eq. (2.46) over that of eq. 
(2.49). These are two ways of interpreting AW that imply different physics. That the 



change in charge can be split as in eq. (2.46) and that the dressed charge is related to the 
dressed energy as W = -f^ is quite convincing. An important factor also is that spin- 
charge separation has been observed experimentally [IS] and to account for it requires an 
understanding of the dressing of charge that extends to non-zero temperatures. Equation 



(2.49) does not provide this. 



2.4 Zero temperature 

Now we turn our attention to the limit of zero temperature. This is a special limit as 



the nature of excitations may change. In section 3J2 the zero temperature limit of the 
Hubbard-Shastry B-model is examined in detail and much of what is indicated here is 
made precise. 

In the limit T — > we see that the functions Y a become singular 



lim T ^ ea(v) < 
lim T ^ € a (v) > <^> 



Hm T ^ Y a (v) 

lim T ^ Y a (v) 





oo 



lim T _> Pa(v) 
lim T ^ p a (v) 



0, 
0. 



(2-51) 



Thus the pseudo-energies e a are better quantities to work with and in terms of them the 



TBA equations (2.11) reduce to 



e Q + e/3 -kQ„ K/3 C 



(2.52) 



Here the subscript on the convolution indicates that the region of integration should be 
restricted to an interval Qp, which is the subinterval of Xp for which ep < 0. We say 
that an a-string is at half-filling if Q a = X a , which implies from (2.51 ) that there are no 



holes in the ground state for such strings. 

First consider the situation when all strings are away from half-filling. Here the 



problem mentioned in the paragraph above eqs. (2.17), (2.18) does not arise and the zero 



temperature limit of the formalism for excitations above equilibrium is straightforward. 



The total change in energy for an excitation, given by (2.30), reduces to 



N 



(2.53) 



as for each string at less than half-filling e™ ax > implies Y^ ax = oo. This is the familiar 
picture in which the pseudo-energies play the role of the dressed energies. 
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Similarly the change in charge (2.46) becomes 



AW = Y, ^ (v aj ) - ) + AW^ ind , (2.54) 

3=1 3=1 

To examine the induced charge, let us denote the boundaries of Q a by ±q a , and so for 
type 1 strings one has Q a = (—q a ,q a ) and 

AWind = ^C a (q a )p a (q a )^a(q a ) - ^rCa(-qa)p a (-q a )^ a (-q a ) ■ (2.55) 

Z7T Z7T 



Recalling (2.19) we see that the induced charge is due to the shift of the rapidities 
corresponding to the boundaries of the intervals Q a . This can be understood as a back- 
reaction of the density, which here at zero temperature occurs at the edge of the Fermi 
sea. For type 2 strings the picture is the same. 

Now we turn to the situation of having some strings in the ground state at half- 
filling, let us say that e™ ax < for some 7-strings. Here one must be careful to only 



consider excitations for which the iV 7 satisfy the selection rules (2.7) as there are no 
holes for 7-strings. Put another way, some of the 7-strings in the ground state may have 
no corresponding mode number in the excited state, due to a decrease in the range of 
mode numbers, and are thus necessarily removed. On the other hand an increase in the 
range of mode numbers will mean that there are some holes in the excited state that do 
not correspond to removed strings and are thus not dynamical. This situation requires 
one to reconsider the nature of the excitations. 

Let us outline a convenient prescription for dealing with excitations that change the 
range of mode numbers of strings which are at half-filling. If the range increases we 
choose to consider only excitations for which all the extra mode numbers are filled. 
In our terminology this means that in such an excitation these extra mode numbers 
always correspond to added strings and thus all holes of the excited state correspond 
to removed strings. Obviously if the range decreases in an excitation then the removed 
mode numbers always correspond to removed strings. We refer to such added and 
removed strings as singular strings. Such singular strings have rapidities that approach 
v max in the limit L — > 00 because they correspond to mode numbers at the edges of 
the range. We refer to the remainder of the added and removed strings as physical 
strings. Note that this prescription does not limit the freedom to capture all possible 
excitations. Indeed any excitation for which not all the extra mode numbers are filled 
can be considered as a limit of an allowed excitation where the rapidities of the necessary 
number of physical removed strings approach t> max . 

An excitation could also result in an overall shift of the mode numbers. This would 
correspond to the removal of some singular strings at one end of the range and the 
addition of singular strings at the other. It can be seen however, that for each of the 
quantities of interest to us, that this transfer of singular strings is not important. 

Let us thus break the added and removed strings into two types, physical and singular 

jV a = j\r 7 pa + jV° a , = iV 7 pr + , (2.56) 



16 



where we use p and s to denote physical and singular respectively. Let us further denote 
the changes in numbers of physical and singular strings as 



5N* = 

Then our prescription is that 



5N* 



(2.57) 



c 7( aoiVi -t- fij 77 '0ivy , (2.58) 

where the right hand side here is the change in the range of mode numbers of 7-strings 
found from eq. (2.6), and we use 7' as a dummy index to make it clear that the sum is 
only over strings which are at half-filling. Let us remark that there may be a restriction 
on the number of physical roots one can excite as only solutions to eq. (2.58) for which 
8N* is an integer for each half-filled string are allowed. 

Now consider again the change in energy formula (2.30) which here takes the form 

7V a 



AE 



3=1 



,1 *| _ pinaxr 



N' 



1 n j 



(2.59) 



Splitting the strings between their physical and singular subsets this becomes 



AE 



N p ' 

E £ ° 



E 



- et 



7 



(2.60) 



where all the constant terms have cancelled due to (2.58). The singular strings may 
also have non-zero dressed momentum and non-trivial dressed scattering. These can 
be redistributed among the physical strings according to the solution of (2.58). For 
example, if the solution to eq. (2.58) is 5N® = f ia 5N^ then the dressed momentum and 
dressed scattering take the following form for the half-filled phase 



>h.f. 



7« > 



(2.61) 



Finally let us consider again the change in charge. As for the energy, the contributions 
of the singular roots cancel all constant terms appearing in the dressed charge (2.47). 
Also the derivative in (2.48) is zero at zero temperature for half-filled strings and so such 
strings do not give rise to an induced charge of the system. The formula for the change 
in charge thus takes the form 



AW 



7V"P a 



N pl 



a.h.f 
ind 1 



where AW^ denotes the induced charge due to the back- reaction of the strings which 
are away from half-filling. 
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2.5 Summary of the excited state formalism 

Let us summarise the main features of the excited state formalism developed in this 
section. In the thermodynamic limit the Bethe equations become the equations for 



densities (2.4). Requiring the free energy to be minimised yields the TBA equations 



(2.11). This closed set of non-linear integral equations on — allows one to determine 

pot 



the equilibrium densities. Excitations of a finite number of strings above the equilibrium 



state can be investigated by examining the shifts of the roots (2.19) arising due to an 
excitation. These shifts also satisfy a closed set of non-linear integral equations (2.20), 
and with the aid of these equations all the features of an excitation can be extracted. 

Let us clear the notations of the previous subsections and parametrise an excited 
state by iV p added particles with rapidities v Pk , k = 1, . . . , N p , and by N h holes with 
rapidities Vh k , k = 1, . . . , N h . Here the indices p^ and h^ include the information of the 
type of string. The excitation can be encoded in a set of Bethe equations for which the 
pseudo-vacuum is the equilibrium state 

ATP ATli 

(-l)v- = e iP Pk( v Pk) L e i *PfcPj(«Pfc> t 'pj) e-^Pfciy^k'"^-) ; (2.62) 
i=i i=i 

N h NP 

(_1)-<p« = e ^K) L Y[ e^^'K'V Y[ e~ i<E WV^ . (2.63) 
Let us stress that these Bethe equations are only valid for large L. The momentum is 



given through eq. (2.36) 



P a = Pa — Pp * (p/3a for particle excitations, (2.64) 
P a = —p a + P/3* (p/3a for hole excitations, (2.65) 



and the scattering phases are determined through the closed set of equations (2.42). 



Due to interactions with the equilibrium state the energy of each excited root gets dressed 



(2.31) 



E a = e a + Tlog(l + ax )fc < g a for particle excitations, (2.66) 
E a = —e a — Tlog(l + max )/c / 3 Q . for hole excitations. (2.67) 
Here e a = T log Y a are the psuedo-energies and the constants and k a $ are defined 



after eq. (2.29). Similarly, for a conserved quantity the corresponding charge of each 



excited root gets dressed and is given by (2.47) 



W a = u a kg a = — for particle excitations, (2.68) 

1 + ymax d^X > K J 



W a = -0J a + , Vmax kp a = — -r-^ for hole excitations, (2.69) 

1 ~r X p O-fJj 
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where u a are the pseudo-charges determined through the closed set of equations (2.44) 
and the constant Co>™ ax is defined above eq. (2.46). The total change in charge for the 
excitation also has a contribution AWind given in eq. (2.48). This is an induced charge 
of the system that is not carried by the excited roots, but rather is due to a back-reaction 
of the densities. 

In the limit of infinite temperature T — > oo the dressed energy and dressed charge 
take their bare values as one expects. The limit of zero temperature T — > requires 
special attention and is discussed in detail in section |2.4| If some strings are at half- 



filling there may be restrictions on allowed excitations and their nature may be altered. 
Regardless of whether this is necessary the dressed energy and dressed charge of excited 
roots take their pseudo values in the zero temperature limit. 



3 Hubbard-Shastry models 



Now we wish to examine the excitations of the Hubbard-Shastry models. An introduc- 
tion to these models can be found in [3], as well as an overview of their equilibrium state. 
We will focus here on the A- and B-models as the excitations of the Hubbard model 
have been investigated elsewhere, see for example [5] and references therein. Not alone 
are these models of great interest but also the B-model in particular provides a good 
illustrative example for the application of the formalism developed in section [2] We will 
restrict our attention to excitations at zero temperature. 

Let us review the Hubbard-Shastry A- and B-models. Here we introduce them al- 
ready in their diagonalised form, their Hamiltonians are presented in detail in [3]. Their 
Bet he equations are [16] 



M 

D iLp(v k ) TT v k 



n 



Vk 



w j + i u 



Wk — Vj — IU t r Wk — W i + ll U 



n 



Wk — Vj + in 



n 



3=1 ■> 1=1 



wi — 2iu 



,N <L, 



M < 



N 



(3.1) 



where e ip ^ for each of the models (and for completeness for the Hubbard model also) are 
given in Table [TJ along with the dispersion relations £ (p). Note that rapidity variable v 
of momentum carrying roots is related to the momentum through y with v = \{y+l/y), 
and so e ip ^ is a double valued function of v. We refer to these roots as ^/-particles and 



they are discussed in detail in appendix A. 3 





A-model 


B-model 


Hubbard 


e ip(v) 

E(p) 


i— 2yu 

-2cosp- 2y/l + u 2 


y+'Mu 
y—2iu 

— 2 cosp 


iy 

— 2cosp — 2u 



Table 1 : Momenta and dispersion relations of Hubbard-Shastry models with v = \ (y + 1/y) . 
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The string hypothesis for the behaviour of the Bethe roots in the thermodynamic 
limit is that each root is a member of one of the following types of strings [TTl Q3] 



• y-particle: a charge 1, spin-up momentum carrying particle, 

• M|uu>-string: a charge 2M, zero spin momentum carrying bound state, 

• M|w-string: a zero charge, spin — M bound state, 

where M denotes a positive integer. The momentum and the dispersion relation for the 
M|vw-strings are 

M 



pm\v W {v) = J2p+( v + & - M H + p-( v - & - M H 

i=i 

M 



(3.2) 



and (Ti\ vw = sign( dp ^"" ) = —1. Here we have adopted the clean notations of [IB] , which 



are compared with the more conventional notations, of say [5], in appendix A.l 



These strings constitute the particle content of the model in the thermodynamic 
limit and they are summarised in figure [T] Re-writing the Bethe equations in terms of 



x : y+ 
x : y- 
• : w 



u{ 



X 
XX 



X 
XX 



XX 
X 



y+ y- 



l\vw 2\vw 3\vw ■ ■ 



l\w 2\w 3\w 



Figure 1: (Colour online) An illustration of string hypothesis for the Hubbard-Shastry models. 
The horizontal line represents the real axis while the imaginary direction extends vertically. 
The y± refer to the two branches of the y-particle and the x mark the corresponding value of the 
t;-rapidity variable, which take values on (—1, 1) for the A-model, and on (— oo, — 1)U (1, oo) for 
the B-model. The circles mark the rapidities of w roots. The M\vw and M\w strings have real 
centres that take values on K. The M|t»w-strings are momentum carrying and their momentum 
Pm\vw an d dispersion £m\vw are obtained by summing the contributions of the roots of which 
they are composed. 
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the string solutions they take the form0 



OO ^M\vw OO ^N\w 

i = e iL p^ n n ^(%* - n n ^(%* - ^ > 

ATy OO ^N\vw 

-l = e^^i-^.^ n ^k* - n n w^,* - ^.o > ( 3 - 3 ) 

3=1 N=l 1=1 

N y OO ^N\w 

1 IT S M {w M ,k - Vyj) TT ] f ? r . 

3=1 N=l 1=1 ' ' 



The S-matrices are given in appendix A.l These are the string Bethe equations which 
represent the starting point, eq. |2.1 of the formalism developed in section [2j 

There are two conserved quantities corresponding to charge and spin, and their re- 
spective chemical potentials are /i and B. The energies carried by the strings are influ- 
enced by these chemical potentials and are 

e y = £y ~ ^ ~ B , &M\vw = £>M\vw ~ 2M/1 , &M\w = 2MB . (3.4) 

We denote the bare charge of an a-string by w° = — and the bare spin by w^ = 



1 ae a 

2 8B 



3.1 A- model 

Now we examine the zero temperature excitations of the A-model. Let us first recall the 
ground state. This is determined by the zero temperature limit of the TBA equations 



(2.11) 



ty £y fl B , 

£m\vw = £m\vw - + e y ® Qy K M , (3.5) 
e M \w = 2MB + t y ® Qy K M ■ 

The dispersions E y and £i\ vw are plotted as functions of v in Figure [2] The y-particles are 
the only strings which have a non-zero density and so the ground state is spin polarised. 
As ^/-particles have charge 1 and spin - the magnetisation of the ground state is half 
its filling. The phase diagram^ is presented in Figure [3] The ground state is empty 
in phase I as the energy cost of having electrons in the ground state is too high. In 
phase II ?/-particles enter the ground state and the filling increases with increasing /x, up 
to the boundary with phase III where the ground state becomes half-filled. Analysing 



the TBA equations (3.5) one can see that the M\v ^-strings are type 1 strings while 



4 These correct a minor error in eqns. (A.11-A.13) of [3J. 

5 We restrict our attention to the quadrant fx < and B > as the other quadrants are related by 
symmetries. 
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Figure 2: A-model: Plots of £ y (v) and £i\ vw (v) for u = 1/2, 1, 2. In the plots of £ y (v) the y+ 
branch is represented by a solid line and the y_ branch is represented by a dashed line. 



Figure 3: A-model: Zero temperature phase diagram in 
the /xl?-plane. The phases identified are: I) empty band, 
II) partially filled and spin polarised band, III) half-filled 
and spin-polarised band. The line separating phases I and 
II is /x + B = 2 — 2\/l + u 2 , and the line separating phases 
II and III is n + B = -2 - 2\/l + u 2 . 




the M|u>-strings are of type 2. The y-particles should be treated as type 1 strings as 
discussed in appendix |A.3| and so the phase shifts are 



<Py,M\vw(v,t) = 7Tby iM \vw + @A/0 ~ *) , <j>y,M\w(v ,t) = 1lb yM \ w + ® M (v ~ t) , 

and the dressed momenta are 

d P y 



P„ 



Py 



p 



M\vw PM\vw 

d P y 



dv 



y 6m — irn y (2by 

,M\vw + 1) 



(3.6) 



M\w 



dv 



y 6m - irn y (2b y! M\w + 1) 



where b y ^u\w and b y ^M\vw are the integers that determine the branch of the excitation. 
Let us first consider excitations above the half-filled phase III. One does not have 



to introduce singular strings through eq. (2.58) as the range of mode numbers for y- 



particles N y + N y = L does not change for an excitation. That k y /3 = for all /3-strings 
is a reflection of this. Moreover the dressed energies can be written explicitly as 



e y ~ £y ~~ ^ ~ P ' e M\vw — — 2M/i , 6m\w — 



M\vw 



+ 2MB. 



while the dressed momenta are branch independent and simplify to 



Py 



M\vw 



0. 



M\w 



-PM\i 



(3.7) 



(3i 



Here the identities (A. 1.8) have been used. The dressed charge and spin are equal to 



their bare values. Thus the y-particles behave as electrons and are gapped, and so there 
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l\vw 



Figure 4: A-model: Plots of e y {P y ) and ex\ vw {P\\ vw ) at u = 1 and B = 0. 

is an energy cost to remove them from the ground state. Their dressed energy and 
momentum take their bare values and so are related as 

e y = -2 cos P y - 2y/l + u 2 - n - B . (3.9) 

The l|w-string is the magnon and the M|w-strings are their bound states. In B = 
magnetic field they have quadratic dispersion at low energies, €m\w ~ ^ 1+ m~ u P~m\ w i 
while in a B > magnetic field they are gapped. The M| ft/;- strings are not dynamical. 
In the strong coupling u — > oo limit the energy gap for removing an electron goes to 
infinity and the physics becomes that of the ferromagnetic spin chair^] [3] . 

Now we consider excitations above phase II where the filling ranges between and 1. 
The dressed energies are given by (3.5 ) and the dressed momentum by ( |3.6[ ), and here the 
choice of branch becomes important. For convenience we choose b y> M\ V w — 0, b y> M\w = — 1 
for all M. The range of dressed momentum is then (— tt, tt) for y-particles, (0, 2it — 2im y ) 
for M\v ^-strings and (0, 2^n y ) for M|w-strings. These should be considered modulo 2n 
and it is convenient to plot the dressed momentum of the M\v if -strings and the M lu- 
strings in the range (0, 2ir). It should be kept in mind that the other branches of the 
excitations are obtained by shifts of 2im y . Plots of dressed energy as functions of dressed 
momentum for ^/-particles, M|w;-strings and M 1 1> unstrings for various filling at B = 
and u = 1 are given in Figures [4] and [5| The velocity can easily be read from the 
plots as it is the derivative of dressed energy with respect to dressed momentum. Away 
from half-filling the dressed charge and spin of some strings gain dependence on rapidity. 
They are given by 

w£ = l, "m\vw = 2M - W M , u c M]w = -W M , 

1 1 1 (3.10) 



6 The ferromagnetic spin-chain Hamiltonian appears at order i and so it is necessary to rescale the 
energies appropiately. 
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Pi 



\w 



Pi 



\w 



Figure 5: A-model: Plots of e 1 | tl) (P 1 | m ), ^i\ w (Pi\ w ) an d ^l\ w ( p i\w) at u = 1 and 5 = 



where Wm (v ) = — 1 ®q K m is a non- negative function that goes to zero both at half- 
filling and zero filling. We see that the ^/-particle behaves as an electron in this phase 
also, but that it is now gapless. In B = magnetic field the M|w-strings retain their 
quadratic dispersion away from half-filling but they are no longer pure spin, they gain 
charge opposite to that of an electron as their energy increases. The dressed spin and 
charge of a 1| unstring is plotted as a function of its dressed momentum for various 
fillings at B = and u = 1 in Figure [5j In a B > magnetic field the M|u>-strings 
become gapped. The M|uu>-strings become dynamical in phase II but they are gapped 
throughout. 



3.2 B-model 

Now we discuss the zero temperature excitations of the B-model. Let us begin again by 
recalling the ground state. The zero temperature TBA equations are 

e y — £y ~ ^ ~ B - €i\ VV] *q 1]vw Ki , 
cm\vw = £m\vw — 2M/J, + e y ®q v Km — £x\ V w *Qi\ vw K\m ? (3-11) 
e M \ w = 2MB + e y ® Qy K M . 



The dispersions £ y and £\\ vw are plotted as functions of rapidity v in Figure |6j Here 
l|vw-strings and y-particles are the only strings that can have non-zero densities. The 
phase diagram is presented in Figure [7j For B = there exist only l|fty-strings and the 
ground state has zero magnetisation. In a B > magnetic field some of the bound pairs 



get broken introducing y-particles to the ground state. From the TBA equations (3.11) 
one can see that the M|vw-strings are type 2 strings while the M|u>-strings are of type 
1, and the ^/-particles should be treated as type 2 strings as discussed in appendix | A. 3 



Let us first focus on excitations over the ground state when it is half-filled and has 
zero magnetisation. This is the subregion of phase V where fio < fi < and B = 0. The 



TBA equations |3.11| can be solved explicitly with solution 

ei\vw — — £ y ®s — fi , e M >2\vw = — 2(M — e y = £ y — £\\ vw *s , e M \ w = 0- (3.12) 



24 



-8 -6-4 -2 2 4 6 8 

V 



-8 -6 -4 -2 2 4 6 8 

V 



Figure 6: Plots of £i\ vw (v) and £ y (v) for the B-model for u = 1/2, 1, 2. In the plots of £ y {v) 
the y+ branch is represented by a solid line and the y_ branch is represented by a dashed line. 



Figure 7: B-model: Zero temperature phase diagram in the 
/U-B-plane. The phases identified are: I) empty band, II) par- 
tially filled and spin polarised band, III) half-filled and spin- 
polarised band, IV) partially filled and partially spin polarised 
band, V) half-filled and partially spin polarised band. The 
value fj,Q ranges between in the limit of weak coupling, and 
—2 + 2 log 2 ~ —0.6137 in the limit of strong coupling. 




The M\v ^-strings for M > 2 and the M|«;-strings are not dynamical. The l|vw-strings 
are half-filled and so we use the prescription (2.58) for handling excitations which change 
the range of mode numbers and get 



5NI 



\vw 



-SNf, 

2 l\vw 



M=2 



P 

M\vw " 



(3.13) 



Thus only excitations with SN^ V w + even are allowed, as SN*, VW must be an integer 



bl\vw,y 



To calculate the dressed momenta let us make the branch choice bi\ VWi M\ V w 
0, 



1 and 



P 'X\vw 
P„ 



-7T +Pl[ t 
-y+Py 



1 ~ dpi|„, 

T * - L - 

27T dv 



Pl\l 



-k S . 



P 



M\w 



~Py ® S 

o, 



P 



M>2\vw 



(M 



where we have introduced the useful functions 
T(i>) = @i * s(v) = i log 
= & 2 *s(v) - 0i (v) 



r( 


i 


r( 


i 




7T 




2 



4u 



ir(i 



'4u> 



2 arctan 



6XP ^ 



l)ir mod 27r , 

(3.14) 

(3.15) 



, , U (|, ^y). Let us remark 



The range of is (— | , |) while the range of P^ is (— y-, 
that when taken modulo 27r the range of P y will have an overlap. The singular l\vw- 
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Pi, p s 

Figure 8: Plots of E h (P h ) and E S (P S ) for the half-filled B-model at n = 0, B = 0. 



strings appearing through eq. (3.13) have rapidity v 
The dressed charge and spin are 



and so carry momentum 0. 



to 



5 U 1\VW ~ 

y 2 ' M\vw 







U M>2\vw 



2M-2 
-M, 



(3.16) 



and thus the phase exhibits spin-charge separation. The removed l|fw-strings get 
dressed as holons while added y-particles get dressed as spinons. The energies and 
momenta of the holons, antiholons and spinons are 



E h 
E s 



-Cl\r 



£y® S + jl. 



i 



Ls<H 



P\x — —Pl\v 
Hi Ph = —Pl\v 



■ Py® S . 
P. 



2\vw 



-Ti sign 



-Py ® S 

-k s . 



(3.17) 



Here an antiholon is identified with composite excitation of a holon and a 2|t>w-string 



because a 2|ww-string is not dynamical. Let us remark that through eqs. (3.12), (3.14) 



and (3.16) it can be seen that its addition can also be regarded as the action of the 
charge su(2) raising operator on a state. Similarly the spinon s is a composite of a 
spinon s and a l|to-string. Plots of Eh(Ph) and E S (P S ) are given in Figure [8] for various 
values of u. The velocities can again easily be read from the slopes. The holons are 
gapped for fi > /i . The gap goes to zero in the weak coupling u — > limit while the 
gap has a maximal value of 2 — 2 log 2 m 0.6137 at \i = in the strong coupling u — > oo 
limit. The spinons are gapless and display an "hourglass" dispersion. The similarity to 
experimental data on spinon scattering in some cuprate materials, say Fig. 2 of [19J, 
cannot go unremarked. In the strong coupling u — > oo limit the two lower wings join 
differentiably at 7r and the upper wings can be understood as the contributions of dressed 
electrons on doubly occupied sites [3]. 

Next let us calculate the dressed scattering phases to examine the scattering of the 
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holons and spinons. First we present explicitly the bare scattering phases 



<f>i\ww,i\vw(v, t) = 2tt - 7T sign(f ) +Q 2 (v-t), 

<I>i\vw,m>2\vw( v > t) = 2tt - 2tt sign(» + Q 1M (v - t) , 

<t>i\vw,y(v, t) = (f>y,i\ vw {v, t) = -ix sign(w) + 9i(v - t) , 

(pi\vw,M\w(v,t) = 0, <f> y ,y{v,t) = 0, </> y ,M\w(v,t) = -7rsign(t;) + Q M (v-t). 

To calculate the dressed scattering it will be necessary to redistribute the contributions 
of the singular strings as in eq. (2.61). Solving eqs. (2.42) for $i\ vw ,i\ vw , &m\vvj, i\vw an d 
® y ,i\vw we get 



$i\vw,i\vw(v, t) = 7i-n sign(w) + T(v) + T(v - t) , 

$M>2\vw,l\vw(v,t) = -27Tsign(w) + Q M -l(v) + &M-l(v ~ t) 

® y ,i\vw(v,t) = -7i - vrsign(w) - - - t) , 



(3.18) 



where the identity tt siga-kKi = ©i has been used. As we are interested in the scattering 
phase shift, which is defined modulo 2n, sign(t>) can be dropped from (3.18). Taking 
into account the scattering phases of singular strings the dressed scattering phases are 



h.f, 

\vw,M>2 



h.f. 



\vw( v >t) =K + @M-l(l> -t) 
,(v,t) = 0, 



^ y ( v ,t) = n-T(v-t) 
^fi\ vw (v,t)=7r + ^(v-t) 

<T) h - f - 

^y,M\w 



h.f. 

y,M>2\vw 



(v,t) =-K + e M (v-t) 

(V,t)=0, 



M\vw,N\w\"i "J "J 

M>2\vw,N>2\vw( V i *) = + &Mn(v ~ t) ~ M +n(v ~ t) ~ M +N-2(v - t) 

(v, t) = 5 MN ir - Q MN (v - t) . 



$ h.f. 



tf,h.f. 

^M\w,N\w 



Let us compute explicitly the phase shifts for the charge triplet and charge singlet 
excitations. 

Charge triplet: holon-holon scattering. Here two l|uw;-strings with rapidities v\ and v 2 

are removed. Let us say that v\ has a greater velocity ^ than v 2 , and let us denote 
this as V\ y v 2 - Then 

and the phase shift is 

5ct = tt + T(«i - v 2 ) ■ (3.19) 

Charge singlet: holon-antiholon scattering. Here two l|vw-strings with rapidities V\ >- v 2 
are removed and a 2|mu-string with rapidity v is added. The rapidity v of the added 
2|t>u>-string can be fixed through eq. (2.41) using P2\ vw = 0, 



2 rw = -$2\vw,l\vv>( V i V l) ~ ®2\v'w,l\vw( V - V 2) + ® 2 \vw,2\vw{ V ~ *0 

= 7T — 6i(f - Vi) - 0i (u - v 2 ) + <d 2 (v - v) , 
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Figure 9: B-model: Plots of dressed energy £i\ vw {P\\ vw ) and dressed charge ^\\ vw {P\\vw) a t 
u = 1 and B = for various fillings. 



and so F2\ vw {v) = n gives v = Vl + V2 . Here 

F llvw {v) = -^iL.ii^^^i) - $?|L,i|^(^,f2) + ®\\L,2\v W (v, Vl+ 2 V2 ) 

and thus 

5cs = T( Vl - u 2 ) - 9x (^p) . (3.20) 

Let us remark that these results as functions of the rapidity are the same as those of the 
half-filled Hubbard model, see e.g. eqs. (7.124) and (7.126) of [5]. The scattering shifts 
for the spin triplet, singlet and spin-charge excitations can be computed similarly and 
also agree with those of the Hubbard model, eqs. (7.139), (7.141) and eq. (7.145) of [5]. 

Now let us consider the less than half-filled phase while still keeping B = 0. In 
Figure [7] this is the portion of phase IV along the /z-axis. The magnetisation is zero and 
the filling is 2n\\ vw . Here again 6m\w = and Pm\w = for M|w-strings but to find the 
dressed energies and momenta of M \v w-strings and ?/-particles one must solve the TBA 
equations (3.11) numerically. Taking all b a/ 3 = 0, the dressed momentum for y-particles 
takes values in irni\ vw < \P y \ < n + irni\ vw , for l|t>u>-strings in irni\ vw < \Pi\ vw \ < vr, for 
M|m/>strings with M > 3 odd in 2Tmi\ vw < \Pm\vw\ < ^ and for M|ww-strings with M 
even in 7r + 2im\\ vw < \Pm\vw\ < 27T. These ranges should be considered modulo 27r but 
it is more convenient in plots to use the ranges specified here. The other branches are 
obtained by shifts of 2nni\ vw . The dressed spins of the excited strings take their bare 
values while the equations for dressed charge are 

ujy = 1 - uj c i1vw *q 1]vv> Ki , u c M \ vw = 2M - u^ vw *q Mvw K 1M , u c M \ w = . (3.21) 

These are rapidity dependent for the ^/-particles and M|fw-strings. Let us remark 
however that at v — ±oo the dressed charges take the values they have at half-filling 



(3.16) 



^(±oo) = 0, ^, w (±oo) = l, ^> 2 | w (±oo) = 2Af-2. (3.22) 
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Figure 10: B-model: Plots of dressed energy e y (P y ) and dressed charge ujy(P y ) at u = 1 and 
B = for various fillings. 



In Figure [9] the dressed energy and dressed charge of a l|t>w-string as function of its 
dressed momentum is plotted for various fillings at u = 1. The corresponding plots 
for y-particles are given in Figure [10] and the "hourglass" behaviour about tt is seen 
again. Both excitations are gapless. The M| fit*- strings with M > 2 are dynamical 
but have a gap of — 2(M — The l|mu-strings remain spinless but y-particles with 
non-zero energy here have dressed charge in addition to their spin. At low energies the 
charge carried by a y-particle scales with its energy, and moreover the magnitude of the 
charge carried increases as the filling is decreased. Thus one does not have spin-charge 
separation away from half-filling. 

It is noteworthy that at B = one can clearly see that the l|t> u>-string is a spin- 
singlet bound state. Let us show this. The spin singlet excitation is achieved by adding 
two y-particles with rapidities V\ and v 2 , and adding a l|w-string with rapidity w which 
we initially take to be arbitrary. The relevant dressed phase shifts are 

®v,v( v >t) = -( K i*Qx\™ ®i\vw, y )(v,t) , 
$i\vw,v(v, t) = -7r sign(w) + 0i (v, t) - (K 2 * Ql|mo $i\ vw , y ) (v, t) , 

®y,l\w(v,t) =TT + Oi(v - t) , ® llw!y (v,t) = 7T + O^V - t) , $l\ Wj l\ w {v , t) = 71 - B 2 (v - t) . 

as P\\ w = 0. The scattering phase 
shift is 

F y (vi) = vi) + ^(vi,^) + vr + 0i(^^) . (3.23) 

Note that 0i(±iu) = ±ioo and so the final term gives rise to a pole of the S-matrix at 
v i — v 2 — 2zti. Although the first two terms cannot be obtained explicitly it can be seen 
analytically that they cancel the pole through the term 0! in $i\ VW)V for v G Qi\ vw , while 
for v ^ Qi\ vw the pole remains. The pole corresponds to a bound state of a y_- and a 
y + -particle as ImP_(f) > for Imv ^ and lmP + (v) < for Imt> ^ 0. The bound 
state, with rapidity v, thus corresponds to a spin singlet excitation with 

Vi = v — iu, v 2 =v + iu, w = v. (3.24) 



The rapidity w is fixed to w = Vl Z" 2 by eq. (2.41 
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Figure 11: B-model: Plots of dressed energy ei\ vw (Pi\ vw ) and dressed charge uj^ vw (P^ vw ) at 
u = 1 and fj, = for various values of magnetisation m. 



and the changes of energy and momentum are 
AE = e llvw (v), AP = P llm 



for v ^ Q 



l\vw 



(3.25) 



Here the identity K\{v + in — iO) + Ki(v — iu + iO) = 5(v) + K 2 (v) was used to obtain 



the energy and eq. (2.36) was used to obtain the momentum. Hence the bound state is 



indeed an added lpty-string. This indicates that spin-spin interactions are responsible 
for the pairing of the electrons into l\v w-strings. 

Finally let us consider briefly the effect of a B > magnetic field at half-filling, 
the interior of phase V in Figure [7j Here there are both l\v u>-strings and y-particles 
in the ground state and it has a magnetisation between and 1/2. The l|vw-strings 
are at half-fillingF] and so excitations must satisfy (3.13), that is, only excitations with 
5Nf, + 5N^ even are allowed. The l|tw-strings are gapped and have dressed charge 1 
while the y-particles have dressed charge zero and are gap less. The M > 2\v ^-strings 
are non-dynamical while the M|u>-strings are dynamical but gapped. The spin dressing 
equations for y-particles and l|vw-strings are 



X\vw 



l\vw 



u s y ® Qy K x - u B 1]vw -k Qllvw K 



1M 



(3.26) 



Let us remark that w^(±oo) = 1, wj, (±oo) = — \ and thus at zero temperature the 
dressed spin jumps as soon as a magnetic field is introduced. This is true for any filling. 



Plots of the dressed energy and dressed spin of a l|mu-string are given in Figure 11 for 
various magnetisations at u = 1 and \x = 0. The corresponding plots for y-particles are 



given in Figure 12 



7 There are no holes for l^ui-strings here. There are less of them than there are at B 
range of mode numbers is decreased by the presence of the y-particles. 



as their 
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Figure 12: B-model: Plots of dressed energy £y{Py) and dressed charge w^(-Py) at u — 1 and 
/i = for various values of magnetisation m. 
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A Appendices 

A.l Conventions, definitions and notations 
Convolutions 

The symbol * denotes the following "convolution" 

POO 

g -k h = / dt g(u,t)h(t,v) , 



(A.1.1) 



If g (or h) is a function of a single variable then one just drops u (or v or both), e.g. if 
g = g(t) then g * h = dt g(t)h(t, v). However, if g or h is a kernel defined through 
a function of one variable then it should be understood as g(t, v) = g(t — v ). 

In section [2] we have expressions of the form g a p * and here convolution is under- 
stood to be over the domain of the rapidity of the /3-string. In section [3] we deal directly 
with the Hubbard-Shastry models for which the ?/-particles have a non-trivial domain 

We use the symbol © to denote a contour integral in 



of rapidity, see appendix A. 3 



the counter-clockwise direction around the branch cut of x(v) } given by eqn. (A. 3. 15). 
Explicitly, for the A-models one has 



g ay © h 



'I*I<1 

while for the B-model 



dt (g a _(u,t)h^(t,v) - g a+ {u,t)h +J (t,v)) = g a „*h^- g a+ *h +1 



g ay © h 



[*I>1 



dt (g a _(u,t)h^(t,v) - g a+ (u,t)h+y(t,v)) = g a - * /i_ 7 - g a + * h +1 
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where 

g a ±(u,t) = g ay (u,t±iO) , h ±1 (t,v) = h yi (t ± iO,v) , 

and * and * denote convolutions with the integration over \t\ < 1 and \t\ > 1 respec- 
tively. 

Matching the notations and conventions 

Most of our notations and conventions come from [18j . and here we compare them to 
those of [5J. 

In the Bethe ansatz we denote particles momenta as pj and auxiliary roots as Wj, so 
they are related to the ones in [5] as pj -H- kj, Wj -B- Aj. 

In the string hypothesis a M|w-string is a A string of length M, a w-particle is a 
A-string of length 1, a M|fw-string is a k-A string of length M, and ?/-particles could 
have been called fc-particles. 

In the TBA equations the Y-functions are related to the ones in [5J as Y M \ W r] M , 
Y M \vw O rj' M , y_(sin(A;)) O ((h) , \k\ < ir/2 and y+(sin(fc)) <B CO) , \k\ > tt/2. 

Kernels and S-matrices 

In section [3] we deal explicitly with various S-matrices and kernels which we list here 
T ~ / \ Id, , , 1 uM , , ?)-iuM 

= 2™ ^ l0g5M(v) = * v^M* ' = ' (A ' L2) 

1 d 



K M n{v) = — -j- logS MN (v) = K m +n{v) + K N - M (v) + 2 V K N _ M+2j {v) , (A. 1.3) 

JW-l 

Smn{v) = S M+N (v)S N ^ M {v) Yl S N -M+2j{v) 2 = S NM (v) , (A. 1.4) 

i=i 

s(u) = ^-log^(u) = - S(u) = -tanh(^-^), (A.1.5) 

w 2m dv y 1 4ucoshf^ w v 4u 4 ; ' v ; 

2u 

Let us give also explicitly the functions Qm( v ) — 2arctan(^), 

M-1 

&mn(v) = &m+n(v) + &n-m(v) + 2 ®N-M+2j(v) . (A.1.6) 

J=l 

Some useful identities 

1 * K M = 1 , if M * K N = K M+N , 1 * s — - , (A 1 7) 

K l - s -k K 2 = s , K M+ i - s* K M — s* K M+2 = . 

d ^®K M = - d ^, £ y ®K M = -£ M]vw . (A.1.8) 
dv dv 1 

7r sign * K M = Q M . (A. 1.9) 
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A. 2 Dressed momentum as a variable 



In the text we have chosen to parametrise all quantities using the f-rapidity variable 
for which the S-matrices take a difference form. Here we would like to point out that 
parametrising in terms of dressed momentum has some interesting features. Dressed 
momentum as a function of v is given in eq. (2.36) as P a = p a — Pp * 4>pa- It is defined 
in terms of the density with respect to which the momentum gets dressed. By inversion 
one obtains v a (P). 

The densities parametrised in terms of P are related to those parametrised in terms 
of v as p a (P) = | | Pa (v a (P)) ■ It is worthwhile here to choose a different normalisation 
for the densities so let us introduce p a = ^Pa and p a = ^p a - The equation for densities 
eq. (2.4) then simplies dramatically to 

p a (P)+p a (P) = l, (A.2.10) 

for the density through which the dressed momentum has been introduced, as can be 
seen from (2.37). Remarkably the densities are non-interacting in this parametrisation. 
Moreover the densities are related to the functions Y a simply as 

1 



Pa{P) 



l + Y a (P) ' 



(A.2.11) 



Next consider the expansion (2.19) in terms of dressed momentum. Let us define £ 
through 

C« (-f«,fc) 



Pa,k Pa,k 

The closed equation for ( a takes the form 



(A.2.12) 



(A.2.13) 



where K a p(P,Q) 



dv a (P) 
dP 



K a p(v a (P) — V/3(Q)), and here the shift function is just — ( a . 



The dressing equations also take a nice form, for example for dressed charge one has 

w a + uppp-kKp a . (A.2.14) 



0J n 



Finally, the induced charge term (2.48) becomes AW ind = — C a oo a * p' a , with p' a 



dP 



A. 3 The ^/-particles of the Hubbard-Shastry models 

In this appendix we describe the y-particles of the Hubbard-Shastry models. These 
do not fit into the classification of strings as type 1 and type 2 in section [2] and so we 
indicate how the formalism developed there can be applied to them. First it is convenient 
to introduce the two functions 



xa(v)=v + v\1 -, xb(v) = v + — v 2 . (A. 3. 15) 
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The function xa has a cut along (—1, 1) and the function x-q has a cut along (—00, — 1) U 
(1, 00). Let us define X+' B and X^' B to be respectively the upper and lower edges these 
cuts. For values of v on the cuts we define xa,b(v) = xa,b(v + *0). Then, since v = 
\{y + 1/y), for any given v one has two possible corresponding y-roots and they can be 
parametrized as 

y+ = x(v), y- = —- } (A.3.16) 
x[v) 

so that in terms of the f-rapidity the |/-roots are split into two subsets. The ^-particles 
can thus be regarded as a string whose f-rapidity takes values on a two-sheeted Riemann 
surface. For the Hubbard and A-models we take x(v) = xa{v) and so for these models 
X y = U X A , while for the B-model we take x{y) = xb(v) and here X y = X+ U I?. 

Now let us address the definition of the counting function for these ?/-particles. As 
for the type 1 and type 2 strings identified in section [2j we wish to define the counting 
function z y so that it is continuous everywhere except at the maximum of the pseudo- 
energy e y of the y-string. This is achieved at the cost of the counting function being 
an increasing function of v, and thus care must be taken about how the densities are 
defined in order to ensure that they are positive. Let us restrict ourselves to the parity 
invariant Hubbard-Shastry models and let us assume without loss of generality] that the 
minimum of e y (v) is on the edge X+' B . Then we define the counting function as 

La + z y (v) = 7T(p y + Lp y (v) + ^ <j) y p (v , vp, n ) , (A.3.17) 

P n=l 

where the range of p is chosen to be continuous everywhere along X^' B except at the 
value v = f max G X_' B , the value of v corresponding to the maximum of e y (v), and the 



scattering phases <f) y p are defined as in eq. (2.13) for the case y + = xa and as in eq. 
(2.14) for the case y + = Xb- For the Hubbard-Shasrty models a + = sign(^ ± ) = — 1 



and cr_ = sign(^j^) = 1 and so z y is an increasing function of v for v E X^' B and it 
is a decreasing function of v for v G X^' B . This ensures that the counting function is 
increasing on the contour clockwise around the cut of xa.,b{v) (that goes along X + ' B on 
the upper side and oppositely along X^' B on the other) with the exception of the point 
^max G X_' at which it is discontinuous. When introducing the densities however it is 
necessary to include the factor cr_ for the density of roots on the edge X^' B in order to 
ensure that this density is positive. Thus we have 



P±+P± 

where 



1 
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dp± 



dv 



+ K ± p* Pf} , (A.3.18) 



K±p = v±K<±($ , fc±p(v) = ^~j~ io Z S ±p( v ) ' (A.3.19) 

8 This is the case for both the A- and B-models. The alternative possibility applies to the Hubbard 
model and proceeds in exactly the same manner under interchange of + and — . 
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and 

Py( v ) = < f s c , Py{v) = < _ ■ (A.3.20) 

\p-[v) it v e x_ Ip-W it t> g x_ 

Finally let us remark that lf +J g = —K_p as S+0 = SL^, and so 

k yl 3 = 1® K y/3 = , (A.3.21) 

for all /3-strings. 
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